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UNDERWATER  ACOUSTIC  BEACON  LOCATION  SYSTEM 

[0001]  The  present  application  claims  the  benefit  of  United 
States  Provisional  Application  Serial  Number  62/297,179  filed  on 
February  19,  2016  by  the  inventor,  Steven  E.  Crocker  and 
entitled  "Underwater  Acoustic  Beacon". 

STATEMENT  OF  GOVERNMENT  INTEREST 
[0002]  The  invention  described  herein  may  be  manufactured  and 
used  by  or  for  the  Government  of  the  United  States  of  America 
for  governmental  purposes  without  the  payment  of  any  royalties 
thereon  or  therefor. 

CROSS  REFERENCE  TO  OTHER  PATENT  APPLICATIONS 
[0003]  None . 

BACKGROUND  OF  THE  INVENTION 
(1)  Field  of  the  Invention 

[0004]  An  acoustic  detection  system  is  provided  in  which  the 
detection  system  can  search  for  and  can  locate  an  underwater 
location  beacon. 
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(2)  Description  of  the  Prior  Art 

[0005]  Various  methods  of  detecting,  locating  and  recovering  an 
Underwater  Location  Beacon  (ULB)  have  been  proposed.  Methods  for 
locating  a  ULB  include  cross-correlating  acoustic  signals  received 
by  two  or  more  hydrophones  in  which  the  hydrophones  are  separated 
at  a  discernable  distance. 

[0006]  In  Raudsep  (United  States  Patent  No.  3,559,161)  a  system 
is  disclosed  where  a  cross-correlation  calculation  is  used  to 
estimate  a  distance  to  an  acoustic  beacon.  While  use  of  the 
system  provides  an  improved  location  estimate  relative  to  that 
provided  by  a  single  hydrophone;  the  system  does  not  provide  a 
signal  gain  or  a  greater  detection  range  than  can  be  achieved  by 
using  an  array  of  spatially  distributed  sensors. 

[0007]  Methods  to  aid  the  search  for  underwater  beacons  such 
aircraft  data  recorders,  include  improvements  to  the  beacons 
themselves.  In  Shope  (United  States  Patent  No.  4,951,263),  a  ULB 
is  disclosed  that  transmits  a  spread  spectrum  acoustic  waveform  as 
opposed  to  time  gated  signals.  While  this  method  can  increase  the 
range  at  which  a  spread  spectrum  signal  transmitted  by  a  ULB  could 
be  detected;  the  proposed  method  would  require  the  replacement  of 
all  underwater  beacons  used  in  all  commercial  aircraft.  This 
would  be  a  considerable  expense. 

[0008]  In  Fabre  and  Bignolais  (United  States  Patent  No. 
8,727,263);  a  method  is  disclosed  where  a  beacon  is  physically 
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relocated  to  a  position  where  the  beacon  is  more  likely  to  be 
detected.  In  particular,  the  beacon  is  expelled  from  an  aircraft 
upon  entry  of  the  downing  aircraft  into  the  water.  The  beacon  is 
designed  to  float  free  of  the  downed  aircraft  and  is  retained  by  a 
lightweight  tether  where  the  beacon  can  ascend  from  the  depth  of 
the  seafloor  to  be  detected  by  a  search  platform. 

[0009]  While  this  method  may  improve  the  probability  of 
locating  a  downed  aircraft;  the  system  needed  for  the  method 
requires  potentially  expensive  modifications  not  only  to  the 
beacon,  but  to  the  downed  aircraft  was  well.  Also,  it  is  unknown 
how  reliable  the  proposed  beacon  deployment  and  tethering  method 
would  be  in  a  real  world  environment. 

[0010]  A  method  to  recover  the  data,  but  not  the  physical 
recorder,  is  disclosed  in  Flynn  (United  States  Patent  No. 
8,005,584).  The  cited  reference  proposes  the  recovery  of  data 
from  a  flight  data  recorder  by  using  an  acoustic  data  link  to 
transmit  data  upon  command  from  a  recovery  platform.  While  the 
method  facilitates  data  recovery  by  eliminating  the  requirement  of 
physical  access  to  the  aircraft  debris  field  on  the  seafloor;  the 
method  does  not  provide  an  advantage  for  detecting  and  locating 
the  aircraft  data  recorder  in  a  first  instance.  The  data  recovery 
method  would  also  require  widespread  and  costly  changes  to 
commercial  airliners. 
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[0011]  While  various  attempts  have  been  made  to  resolve  the 
problem  of  detecting  and  locating  a  ULB  on  the  seafloor;  the 
proposed  improvements  require  costly  and  widespread  modifications. 
As  such,  there  is  a  continuing  need  for  a  system  and  method  of  use 
that  increases  the  distance  at  which  a  ULB  can  be  detected;  yields 
a  precise  estimate  of  the  beacon  location;  and  does  not  require 
costly  modifications  to  the  aircraft,  data  recorders  or  acoustic 
beacons  attached  to  the  recorders . 

[0012]  Underwater  Location  Beacons  (ULB)  of  the  type  to  which 
this  invention  pertains,  are  typically  affixed  to  cockpit  voice 
and  flight  data  recorders  used  in  commercial  aircraft.  The  ULB 
presently  used  is  a  battery  powered  transducer  that  transmits  an 
acoustic  signal.  The  acoustic  signal  is  intended  to  be  received 
and  to  aid  a  search  and  recovery  effort  following  the  entry  of  an 
aircraft  into  a  body  of  water. 

[0013]  When  an  aircraft  equipped  with  a  data  recorder  enters 
the  water;  the  ULB  is  energized  by  a  water  activated  battery. 

Once  energized,  the  ULB  transmits  an  acoustic  waveform  at  a  fixed 
interval  of  approximately  one  second  until  the  battery  is 
exhausted  and  the  beacon  ceases  to  operate.  Since  a  precise 
estimate  of  the  location  where  an  aircraft  enters  the  water  may 
not  be  available;  the  search  area  can  be  large  with  the  area 
sometimes  measured  in  thousands  of  square  kilometers.  Due  to  the 
limited  operating  life  of  the  ULB  and  the  potential  for  a  large 
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search  area;  a  ULB  detection  system  preferably  has  the  following 
features . 

[0014]  The  ULB  detection  system  must  be  able  to  be  rapidly 
transported  to  the  location  where  the  aircraft  was  estimated  to 
have  entered  the  water  and  must  be  rapidly  placed  into  operation 
so  as  to  maximize  the  time  to  search  for  the  ULB  prior  to 
depletion  of  the  ULB  battery.  Also,  the  system  must  maximize  the 
rate  at  which  an  area  of  the  seafloor  can  be  searched  for  signals 
transmitted  by  a  ULB  and  the  system  must  be  able  to  detect  the  ULB 
signal  at  the  greatest  practical  distance. 

[0015]  Furthermore,  the  system  must  produce  an  estimate  of  the 
ULB  location  to  aid  in  recovery  of  the  ULB  and  the  flight  data 
recorder  to  which  the  ULB  is  attached.  Also,  the  system  must  be 
capable  of  detecting  and  estimating  the  position  of  the  ULB 
located  in  water  depths  ranging  from  arbitrarily  shallow  to  full 
ocean  depth  in  the  region  of  an  abyssal  plain. 


SUMMARY  OF  THE  INVENTION 

[0016]  Accordingly,  it  is  a  primary  object  and  general  purpose 
of  the  present  invention  to  provide  a  system  to  search  for  an 
Underwater  Location  Beacon  (ULB)  in  varying  water  depths. 

[0017]  It  is  a  further  purpose  of  the  present  invention  to 
provide  a  method  of  use  for  a  linear  acoustic  sensor  array  in 
which  the  method  produces  a  precise  estimate  for  the  location  of  a 
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ULB  using  data  collected  during  a  single,  non-maneuvering  pass  of 
a  typical  search  pattern. 

[0018]  It  is  a  still  further  object  of  the  present  invention  to 
detect  the  presence  of  a  signal  transmitted  by  the  ULB  and  to 
provide  an  unambiguous  estimate  of  the  ULB  on  the  seafloor  from 
which  that  signal  originated  during  a  single  pass  of  the  tow  body. 
[0019]  To  attain  the  objects  of  the  present  invention,  a  system 
and  method  of  use  is  disclosed  that  includes  an  underwater  body 
that  can  be  towed  behind  a  ship.  The  underwater  tow  body  is 
equipped  with  a  linear  sensor  array  comprising  a  plurality  of 
hydrophones.  The  tow  body  includes  electronics  to  provide 
electric  power  to  the  sensor  array  and  components  to  transmit 
electrical  signals  from  the  sensor  array  to  a  data  processor 
located  on  the  ship.  The  processor  uses  the  data  provided  by  the 
[0020]  tow  body  to  detect  and  estimate  the  position  of  an 
acoustic  signal  transmitted  by  the  ULB. 


BRIEF  DESCRIPTION  OF  THE  DRAWINGS 
[0021]  Other  objects,  features  and  advantages  of  the  present 
invention  will  be  apparent  with  reference  to  accompanying  drawings 
in  which  is  shown  an  illustrative  embodiment  of  the  invention, 
wherein  corresponding  reference  characters  indicate  corresponding 
parts  throughout  the  several  views  of  the  drawings  and  wherein: 
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[0022]  FIG.  1A  is  a  perspective  view  of  an  underwater  tow  body 
of  the  present  invention; 

[0023]  FIG.  IB  is  a  top  view  of  the  underwater  tow  body  of  the 
present  invention; 

[0024]  FIG.  1C  is  a  front  view  of  the  underwater  tow  body  of 
the  present  invention; 

[0025]  FIG.  ID  is  a  side  view  of  the  underwater  tow  body  of  the 
present  invention; 

[0026]  FIG.  2  is  a  diagram  identifying  major  components  of  the 
system  of  the  present  invention; 

[0027]  FIG.  3  is  a  depiction  of  information  collected  by  the 
system  as  displayed  to  an  operator; 

[0028]  FIG.  4  is  a  geometry  of  a  linear  sensor  array  of  the 
present  invention  and  a  system  of  coordinates  associated  with  the 
sensor  array; 

[0029]  FIG.  5  depicts  the  tow  body  and  linear  sensor  array  of 
the  present  invention  with  spatial  response  patterns  of  the  array 
shown; 

[0030]  FIG.  6  is  an  illustration  of  the  method  of  use  of  the 
present  invention  wherein  a  tow  body  searches  the  seafloor  for  an 
underwater  location  beacon; 

[0031]  FIG.  7  is  an  illustration  of  the  tow  body  and  the 
spatial  coordinate  systems  used  to  estimate  the  position  of  the 
underwater  location  beacon  on  the  seafloor;  and 
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[0032]  FIG.  8  is  an  illustration  of  the  position  estimate  for 
the  underwater  location  beacon  located  on  the  seafloor  with  the 
position  presented  in  geographic  coordinates. 


DETAILED  DESCRIPTION  OF  THE  INVENTION 
[0033]  A  known  passive  sonar  equation  is  used  to  estimate  the 
maximum  distance  at  which  an  underwater  tow  body  is  likely  to 
detect  an  acoustic  signal  transmitted  by  the  Underwater  Location 
Beacon  (ULB) .  The  Figure  of  Merit  (FOM)  equation  for  a  passive 
sonar  detection  system  is  defined  as  the  maximum  attenuation  in 
the  acoustic  intensity  of  the  transmitted  signal  that  can  occur 
for  successful  detection  of  the  transmitted  signal.  The  FOM  is 
determined  in  decibels  as 


[0034] 


FOM  =  SL-NL+DI -DT 


(i) 


where  "SL"  is  the  source  level  of  the  ULB  signal  in  decibels 
referenced  to  one  micro-pascal  at  one  meter  from  the  acoustic 
center  of  an  acoustic  transducer  (dB  re  1  uPa@m) ;  "NL"  is  the 
spectrum  level  of  interfering  noise  in  decibels  referenced  to 
one  micro-pascal  squared  in  a  one  hertz  band  (dB  re  1  uPa2/Hz)  , 
"DI"  is  the  directivity  index  of  the  acoustic  sensor  in  decibels 
and  " DT "  is  the  detection  threshold  of  the  acoustic  sensor  in 
decibels . 


8  of  31 


Attorney  Docket  No.  300087 

[0035]  Sonar  performance  parameters  that  can  apply  in  a  search 
for  an  aircraft  data  recorder  ULB  in  the  deep  ocean  are:  SL  =  160 
dB  re  1  uPa@lm;  Wi  =  45  dB  re  1  uPa2/Hz;  DI  =  3  dB;  DT  =  25  dB; 
and  FOM  =  88  dB .  The  provided  parameters  invoke  a  number  of 
simplifying  assumptions  including  moderate  weather  conditions  with 
wind  speeds  of  five  to  ten  meters  per  second,  an  acoustic  sensor 
with  a  hemispherical  directional  response  and  a  simple  square  law 
detector  providing  a  conservative  estimate  for  the  detection 
threshold  of  a  typical  ULB  signal. 

[0036]  The  attenuation  of  the  ULB  signal  results  from  a 
geometric  divergence  which  is  assumed  to  be  the  result  of 
spherical  spreading.  Absorption  of  sound  by  seawater  is  also  an 
important  factor  at  the  ULB  operating  frequency  of  37.5  kHz,  with 
an  absorption  coefficient  of  approximately  a  =  0.005  decibels  per 
meter  (dB/m)  near  the  seafloor  in  the  deep  ocean.  Thus,  the 
transmission  loss  TL  is  computed  as 

[0037]  TL  =  201oq10R+aR  (2) 

where  " R "  is  a  distance  from  an  underwater  location  beacon  to  a 
linear  sensor  array. 

[0038]  The  maximum  distance  of  a  known  detection  of  a  ULB  is 
3.5  kilometers.  This  distance  is  the  predicted  detection  range  R 
based  on  consideration  of  Equations  (1)  and  (2)  together  with  the 
parameters  for  SL,  NL,  DI,  DT  and  alpha.  Any  sensor  system  for 
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detecting  a  ULB  signal  must  operate  in  relatively  close  proximity 
to  the  seafloor  in  order  to  minimize  signal  attenuation  and  to 
increase  the  likelihood  of  detecting  the  acoustic  signal 
transmitted  by  the  ULB. 

[0039]  Referring  now  to  FIG.  1A  thru  FIG.  ID,  the  present 
detection  system  comprises  a  tow  body  10  attachable  to  a  tow  cable 
100  by  a  pad  eye  12  on  a  top  surface  of  the  tow  body.  The  tow 
cable  100  provides  an  attachment  to  a  surface  ship  (not  shown) 
used  to  pull  the  tow  body  10  through  the  water.  In  addition,  the 
tow  cable  100  includes  electrical  conductors  that  provide 
electrical  power  or  a  coaxial  cable  (not  shown)  to  support  a 
digital  data  link  between  the  surface  ship  and  the  tow  body  10. 

The  tow  body  10  also  comprises  a  vertical  stabilizer  14  for  stable 
orientation  of  the  tow  body  when  towed  through  the  water. 

[0040]  The  underside  of  the  tow  body  10  supports  a  bracket  16 
for  mechanical  attachment  of  a  pressure  vessel  18  and  an  acoustic 
sensor  array  20  to  the  tow  body.  The  bracket  16  is  shaped  to 
provide  a  set  of  skids  17  that  support  the  tow  body  10  when  placed 
on  a  horizontal  surface,  such  as  the  deck  of  a  ship,  and  to 
prevent  damage  to  the  sensor  array  20  when  the  sensor  array  is 
attached  to  the  tow  body. 

[0041]  The  acoustic  sensor  array  20  comprises  a  linear 
arrangement  of  hydrophones  (not  shown)  oriented  perpendicular  to  a 
longitudinal  axis  30  of  the  tow  body  10.  Thus  the  hydrophones  are 
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generally  orthogonal  to  the  direction  in  which  the  tow  body  moves 
through  the  water  when  towed.  In  the  present  invention, 
orthogonal  is  defined  as  the  longitudinal  axis  of  the  acoustic 
sensor  array  20  being  perpendicular  to  the  movement  of  the 
direction  of  movement  of  the  tow  body  10.  The  tow  body  10  is 
optionally  fitted  with  an  underwater  location  beacon  (ULB)  40  that 
can  be  activated  if  the  tow  body  becomes  detached  from  the  tow 
cable  100  and  when  the  ULB  is  used  (if  needed)  to  aid  in  recovery 
of  the  tow  body. 

[0042]  The  tow  body  10  also  includes  a  magnetic  heading  sensor 
41  to  determine  the  direction  of  the  longitudinal  axis  30  of  the 
tow  body  with  respect  to  magnetic  north.  The  tow  body  10  further 
includes  a  gravity  sensor  22  to  measure  an  angular  orientation  of 
the  tow  body  relative  to  the  gravity  of  the  Earth  and  a 
hydrostatic  pressure  sensor  (collocated  with  the  heading  sensor 
41)  to  indicate  depth  beneath  the  ocean  surface. 

[0043]  Referring  to  FIG.  2,  data  collected  by  sensors  of  the 
tow  body  10  are  received  by  an  electronic  assembly  contained  in 
the  tow  body  pressure  vessel  18.  Electronic  circuits  of  the 
assembly  receive  and  condition  the  analog  sensor  signals  prior  to 
a  digital  conversion.  Digital  sensor  signals  are  formatted  for 
transmission  to  the  data  processor  using  well-known  digital 
communication  standards. 


11  of  31 


Attorney  Docket  No.  300087 

[0044]  A  telemetry  interface  transmitter  of  the  tow  body  10 
sends  the  digital  data  over  the  tow  cable  100  where  the  cable  is 
passed  from  a  deployment  winch  102  to  a  data  processor  200  on  the 
surface  ship.  The  data  is  then  received  by  a  telemetry  interface 
receiver  and  stored  to  memory  to  be  processed  by  acoustic  and 
geospatial  processing  algorithms  (detailed  further  in  the  present 
application) .  Data  processing  results  are  shown  on  an  operator 
display  202. 

[0045]  As  shown  in  FIG.  3,  the  operator  display  202  provides 
information  to  aid  in  the  search  effort  and  to  supervise  movement 
of  the  tow  body  10.  The  operator  display  202  includes  a 
geographic  summary  of  the  search  and  ULB  detection  history  in  a 
geodetic  referenced  display  204.  The  attitude  and  status  of  the 
tow  body  10  includes  roll,  pitch,  yaw,  depth  below  surface  and 
altitude  above  the  seafloor  as  shown  in  display  206.  The  acoustic 
power  spectrum  observed  by  a  user  selected  sensor  array  beam  is 
shown  on  display  208.  The  acoustic  power  observed  in  a  user 
specified  frequency  band  is  presented  to  the  operator  as  a 
function  of  time  and  beam  angle  in  a  three  dimensional  "waterfall" 
display  210. 

[0046]  As  shown  in  FIG.  4,  the  present  invention  uses  the 
linear  sensor  array  20  of  acoustic  pressure  sensors  (hydrophones) 
22  arranged  in  alignment.  This  arrangement  improves  the  maximum 
distance  at  which  the  signal  transmitted  by  a  ULB  is  likely  to  be 
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detected.  The  linear  sensor  array  20  provides  an  improved 
directivity  index  DI  for  improving  the  Figure  of  Merit  ( FOM )  of 
the  detection  system  -  See  Equation  (1) .  The  directivity  index  is 
improved  by  the  spatial  response  of  an  array  beam  pattern  -  as 
further  described  by  the  use  of  Equation  (3)  and  Equation  (5) . 
[0047]  A  Cartesian  system  of  coordinates  is  defined  such  that 
the  pressure  sensors  (hydrophones)  22  are  arranged  in  a  line  along 
the  "Ya"  axis  with  "d"  as  the  spacing  between  adjacent  hydrophones 
(subscript  "a"  indicates  that  the  system  of  coordinates  is 
referenced  to  the  linear  sensor  array  20) .  Thus,  the  Ya 
coordinate  points  in  the  starboard  direction  are  also  relative  to 
the  tow  body  10.  The  Xa  coordinate  points  forward  and  the  Za 
coordinate  points  in  the  direction  defined  by  a  right  hand 
coordinate  system.  In  straight  and  level  flight,  the  Za 
coordinate  would  point  to  the  seafloor. 

[0048]  In  the  figure,  a  system  of  spherical  coordinates  is 
defined  where  the  coordinates  of  an  arbitrary  point  Pa  in  the 

field  is  located  at  an  azimuth  angle  in  the  Xa-Ya  plane  of  yfa  ,  a 

zenith  angle  from  the  Za  axis  of  ,  and  at  a  distance  r  from  the 
origin  (the  acoustic  center  of  the  linear  array  20) .  The  range  of 
y/a  is  +/- 90°,  while  the  range  of  is  +/-1800.  The  angle  is 
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measured  in  the  Xa~Ya  plane  with  the  reference  line  of  measurement 


being  the  Xa  axis  while  "a  is  measured  from  the  Za  axis  to  the 


line  (vector)  that  points  from  the  acoustic  center  of  the  linear 
array  20  to  the  point  Pa  in  the  field. 

[0049]  By  way  of  example;  the  linear  sensor  array  20  typically 
comprises  thirty-two  individual  hydrophones  22  spaced  at  an 
interval  of  fifteen  millimeters;  thereby,  resulting  in  an  aperture 
length  of  approximately  46.5  centimeters.  A  beam  pattern  b  for 
the  described  linear  sensor  array  300  is  calculated  using  well- 
known  methods  [L.J.  Ziomek,  Underwater  Acoustics :  A  Linear  Systems 
Approach ,  Academic  Press  (1985) ]  as 


where  N  is  the  number  of  hydrophones  22  in  the  linear  sensor 
array  300;  n  is  an  index  used  to  identify  an  individual 
hydrophone  and  as  the  index  of  summation  for  Equation  (3)  where 
the  index  varies  from  1  to  N/2;  d  is  the  distance  between 
adjacent  hydrophones,  and  w  is  an  amplitude  weighting  function 
applied  to  the  hydrophones  .  In  Equation  (3)  ,  ky  is  the 
component  of  the  acoustic  wavenumber  vector  that  is  parallel  to 
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the  longitudinal  axis  of  the  linear  sensor  array  20  and  k  "hat" 
is  the  component  of  the  acoustic  wavenumber  vector  that  is 
parallel  to  the  longitudinal  axis  of  the  linear  sensor  array  to 
which  a  given  beam  is  steered.  These  values  are  related  to 
angles  in  spherical  coordinates  in  accordance  with  Equation  (4) . 
[0051]  The  present  analysis  employs  a  Taylor  window  function 
[E.  Brookner,  Practical  Phased  Array  Antenna  Systems ,  Artech  House 
(1991)],  but  other  amplitude  windows  are  not  precluded  from  use. 


The  response  of  the  linear  sensor  array  20  in  the  direction 


a>  a 


is  dependent  on  the  direction  to  which  the  main  response  axis  of 


the  beam  pattern  b  is  steered  ¥a>“a  such  that 


[0052] 


i  .  Q  .  *  2tz  ,  ;  ,  „ 

ky~~j~  &a  ¥a  and  ^ySin0asin^ 


(4) 


with  A  being  the  acoustic  wavelength. 

[0053]  Equation  (4)  is  used  to  calculate  the  beam  patterns 
using  Equation  (3) .  Given  the  beam  patterns,  the  Directivity 
Index  DI  is  calculated  as  needed  to  predict  the  detection  range  R 
and  to  show  the  usefulness  of  the  method  of  the  present  invention. 
[0054]  A  computationally  efficient  method  for  beamforming  the 
linear  sensor  array  20  is  based  on  well-known  frequency-wavenumber 
methods  as  taught  by  the  prior  art  [J.  Capon,  "High-Resolution 
Frequency-Wavenumber  Spectrum  Analysis",  Proceedings  of  the  IEEE, 
vol .  57  (1969)].  The  frequency-wavenumber  beamformer  transforms  N 
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data  channels  of  time  series  data  of  M  samples  each  into  N 
frequency  spectra.  Each  hydrophone  channel  produces  a  time  series 
of  voltage  versus  time.  The  Fourier  Transform  of  the  time  series 
data  is  taken  such  there  are  N  frequency  spectra. 

[0055]  For  this  invention,  a  frequency  spectrum  is  defined  as  a 
representation  of  the  acoustic  data  as  a  function  of  frequency 
versus  time.  As  stated,  the  frequency  spectra  is  computed  using 
the  Fourier  transform  of  the  time  series  data  provided  by  a 
hydrophone.  The  hydrophone  represents  the  signal  as  a  voltage 
spectrum  (volts  squared  versus  frequency) . 

Each  N  data  channel  results  in  one  frequency  spectrum;  thus 
there  are  N  spectra.  M  is  the  number  of  data  samples  in  each 
data  channel.  Each  frequency  spectrum  is  associated  with  one 
beam  steering  angle,  wherein  the  steering  angles  are  uniformly 

spaced  in  the  cosine  of  the  azimuth  angle  yfa  . 

[0056]  FIG.  5  illustrates  beam  patterns  400  formed  by  a 
frequency-wavenumber  beamformer  when  applied  to  data  generated  by 
the  linear  sensor  array  20.  A  well-known  property  of  linear 
sensor  arrays  is  that  the  beam  patterns  b  are  axially  symmetric 
about  a  longitudinal  line  that  passes  through  the  individual 
sensors.  This  arrangement  results  in  a  conical  ambiguity  for  an 
arrival  direction  of  an  acoustic  signal  received  by  the  linear 
sensor  array  20. 
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[0057]  In  a  known  arrangement  where  the  axis  of  the  linear 
sensor  array  20  is  co-linear  with  a  direction  of  movement;  the 
effect  is  often  referred  to  as  a  "left-right  ambiguity"  since  it 
is  generally  assumed  that  the  signal  arrived  from  a  location  in 
the  horizontal  plane,  but  it  is  not  possible  to  resolve  whether 
the  signal  arrived  from  the  left  or  the  right  relative  to  a  course 
heading  of  the  deploying  surface  ship. 

[0058]  Note  that  the  linear  sensor  array  20  employed  by  the 
present  invention  is  oriented  perpendicular  to  the  direction  of 
movement  of  the  tow  body  10.  Therefore,  the  conical  beam  patterns 
are  symmetric  about  the  Ya  coordinate;  oriented  in  the  port  and 
starboard  direction.  The  present  invention  results  in  a  forward- 
aft  ambiguity  where  the  side  (port  or  starboard)  from  which  a 
signal  arrives  is  known  unambiguously,  but  the  forward  or  aft 
direction  is  not. 

[0059]  This  uncertainly  is  illustrated  in  the  figure  where 
select  beam  steering  index  values  are  labeled.  For  example,  the 

beam  labeled  ¥a\+ 1]  is  a  broad  cone  that  intersects  the  horizontal 
plane  at  a  few  degrees  to  starboard  in  the  forward  and  aft 
directions.  Likewise,  beam  index  Va  [+8]  is  further  to  starboard, 
but  is  ambiguous  in  the  forward  and  aft  directions.  Beam  index 

[-12  ]  is  oriented  to  port,  and  is  also  ambiguous  in  the  forward 
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and  aft  directions.  Only  beams  Wa  [±16  ]  are  unambiguous  and  are 

commonly  referred  to  as  "end-fire  beams"  in  that  the  beams  are 
oriented  along  an  axis  of  symmetry  for  the  linear  sensor  array  20. 
It  is  noted  that  FIG.  5  uses  a  view  beneath  the  tow  body  10  such 
that  an  observer  would  be  looking  up  at  the  tow  body.  The  beams 
with  a  positive  index  are  to  starboard  and  the  beams  with  a 
negative  index  are  oriented  to  port. 

[0060]  The  maximum  distance  at  which  a  ULB  signal  can  be 
detected  is  increased  by  the  use  of  the  linear  sensor  array  20. 

In  particular,  the  Figure  of  Merit  ( FOM )  for  the  detection  system 
is  increased  as  a  result  of  the  increased  directivity  index  DI  of 
the  linear  sensor  array  20,  as  shown  by  the  use  of  Equation  (1) . 

As  is  known  in  the  art,  the  directivity  index  is  increased  by  the 
directionality  of  the  linear  sensor  array  20.  The  directivity 
index  of  any  array  is  computed  by  using  well-known  Equation  (5) . 
[0061]  In  Equation  (5) ,  the  directivity  index  DI  for  the  beam 

pattern  developed  by  an  array  of  sensors  is  defined  as  the  ratio 
of  the  power  output  by  an  omnidirectional  sensor  located  in  an 
isotropic,  incoherent  noise  field  to  the  beamformed  output  of  the 
linear  sensor  array  20  in  the  same  noise  field,  given  as 

f  \ 


[0062] 


DI  =  10  log10 


(5) 
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where  dQ  is  the  differential  element  of  a  solid  angle.  The 
differential  element  is  the  differential  used  to  define  the 
integral.  The  difference  element  represents  a  solid  angle  in 
three-dimensional  space.  There  are  4  pi  radians  in  a  sphere. 
[0063]  The  directivity  index  DI  for  the  linear  sensor  array  20 
considered  by  the  present  invention  is  approximately  13  dB  at  37.5 
kHz  (frequency  of  operation  for  a  typical  ULB) .  Therefore,  under 
the  same  set  of  conditions  that  were  employed  to  estimate  the 
Figure  of  Merit  (FOM)  for  a  known  beacon  location  system,  a 
similar  system  using  the  thirty-two  element  linear  sensor  array  20 
would  provide  a  FOM  of  98  dB,  a  10  dB  improvement. 

[0064]  Beam  patterns  b  are  computed  using  Equation  (3) .  The 
directivity  index  for  each  beam  is  computed  using  Equation  (5) , 
which  in  the  provided  example  is  approximately  13  dB .  The 
directivity  index  of  an  existing  beacon  is  estimated  to  be  3  dB 
based  on  the  stated  "hemispherical  response"  of  the  hydrophone 
where  noise  arriving  from  above  is  assumed  to  be  shielded  by  the 
tow  body  10.  Thus,  the  improvement  is  lOdB  (the  difference 
between  the  two) .  The  resulting  maximum  detection  range  for  the 
system  described  by  this  invention  is  4.8  kilometers,  a  1.3 
kilometer  improvement  over  standard  practices  in  the  same 
operating  conditions.  For  example:  current  pinger  locators  do  not 
use  linear  arrays. 


19  of  31 


Attorney  Docket  No.  300087 

[0065]  Another  useful  comparison  between  present  detection 
practices  and  the  system  described  herein  is  the  total  area  of  the 
seafloor  from  which  a  ULB  signal  is  likely  to  be  detected. 

Assuming  that  the  tow  body  10  is  operated  five  hundred  meters 
above  a  level  seafloor;  the  area  beneath  the  tow  body  for  which  a 
ULB  signal  is  likely  to  be  detected  in  present  practices  is 
thirty-eight  square  kilometers.  A  search  system  equipped  with  the 
linear  sensor  array  20  would  be  equally  likely  to  detect  a  ULB 
over  a  seafloor  area  totaling  seventy-two  square  kilometers  - 
nearly  twice  the  area  covered  by  present  beacon  location  systems. 
[0066]  Data  provided  by  the  linear  sensor  array  20  is  also 
processed  to  provide  a  precise  estimate  of  the  specific  location 
of  a  detected  ULB.  As  a  result  of  the  symmetry  in  beam  patterns 
formed  by  the  linear  sensor  array  20;  the  most  likely  location  for 
the  origin  of  a  detected  ULB  signal  is  defined  by  the  intersection 
of  the  conical  pattern  of  the  beam  (in  which  the  signal  was 
detected)  with  the  seafloor. 

[0067]  Thus,  knowledge  of  the  location  and  altitude  of  the  tow 
body  10,  combined  with  the  beam  steering  angle  in  which  a  ULB 
signal  was  detected,  can  be  used  to  identify  the  region  on  the 
seafloor  where  the  ULB  is  most  likely  to  be  found.  The 
intersection  of  the  conical  beam  with  the  seafloor  is  defined  by  a 
conic  section,  in  particular  a  hyperbola.  The  area  about  the 
hyperbola  that  may  contain  the  ULB  is  defined  by  the  width  of  the 
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main  response  axis  for  the  various  beam  steering  angles  generated 
by  the  beamformer.  As  such,  the  ULB  can  be  located  by  many  beams 
as  the  beams  in  which  the  ULB  is  located  will  change  as  the  tow 
body  passes  by  the  ULB  and  the  angle  from  the  longitudinal  axis  of 
the  linear  array  20  to  the  ULB  changes. 

[0068]  FIG.  6  illustrates  the  projection  of  the  conical  beam  in 
which  a  ULB  signal  was  detected  onto  the  seafloor.  The  total  area 
of  seafloor  500  on  which  a  ULB  signal  is  likely  to  be  detected  is 
illustrated  by  the  circular  area  beneath  the  tow  body  10. 

However,  the  likely  location  of  the  ULB  (and  associated  aircraft 
wreckage  600)  is  confined  to  the  subset  of  the  total  surveillance 
area  defined  by  intersection  A  of  the  conical  beam  in  which  the 
ULB  signal  was  detected  with  the  seafloor  500. 

[0069]  As  shown  in  the  figure,  the  directional  response  of  the 
linear  sensor  array  20  identifies  the  ULB  location  on  the 
starboard  side  of  the  tow  body  10.  However,  the  ULB  (and  aircraft 
wreckage  600)  could  be  located  anywhere  within  the  hyperbolic 
projection  of  the  detection  beam  onto  the  seafloor.  Based  on  a 
single  detection,  the  ULB  is  equally  likely  to  be  located  ahead 
of,  on  the  starboard  beam,  or  behind  the  tow  body  10.  However, 
when  the  geolocation  data  from  multiple  ULB  detections  acquired 
during  a  single  non-maneuvering  pass  of  the  tow  body  10  are 
accumulated;  a  precise  location  estimate  of  the  ULB  location  can 
be  developed.  The  accumulation  of  detections  is  based  on  the  ULB 
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"pinging"  at  a  rate  of  once  per  second  -  every  ping  represents  a 
detection  opportunity. 

[0070]  The  area  of  the  seafloor  500  on  which  a  detected  ULB  is 
likely  to  be  located  is  generated  by  projecting  the  conical  beam 
pattern  defined  in  a  moving  spherical  coordinate  system  that  is 
relative  to  the  tow  body  10  onto  a  stationary  Cartesian  coordinate 
system  defined  relative  to  the  Earth. 

[0071]  FIG.  7  illustrates  the  coordinate  systems  employed  by 
the  present  invention.  A  local  North-East-Down  (NED)  coordinate 
system,  as  commonly  used  to  define  navigation  and  ground 
coordinates  is  defined  relative  to  an  arbitrary  point  on  an  ocean 
surface  used  as  a  reference  for  the  search  area.  The  system  is 
defined  by  three  cardinal  directions  as  geodetic  North  (xn) , 
geodetic  East  (yn)  and  Down  (zn)  relative  to  the  World  Geodetic 
System  (WGS  84)  ellipsoid  model  [NIMA  Technical  Report  TR8350.2, 
Chapter  3  "Department  of  Defense  World  Geodetic  System  1984,  Its 
Definition  and  Relationships  with  Local  Geodetic  Systems"  (2000) ] . 
[0072]  Coordinate  vectors  expressed  in  the  local  NED  coordinate 
system  are  denoted  with  the  subscript  "n" .  The  position  vector  Pn 
for  a  location  in  the  local  NED  system  is  denoted  by 


[0073] 


yn 


(6) 
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[0074]  A  vehicle  carried  NED  coordinate  system  is  associated 
with  the  tow  body  10.  The  origin  is  defined  as  the  acoustic 
center  of  the  linear  sensor  array  20.  The  vehicle  carried  NED 
coordinate  directions,  denoted  by  the  subscript  "nv" ,  are  parallel 
to  the  local  NED  coordinates.  Thus,  the  x-axis  (xnv)  points 


toward  geodetic  North,  the  y-axis  (ynv)  points  East,  and  the  z- 


axis  ( znv)  points  down.  The  position  vector  Pnv  for  a  location  in 


the  vehicle  carried  NED  system  is  denoted  by 


[0075] 


Xnv 

ynv 

Znv 


(7) 


[0076]  The  array  (or  acoustic)  coordinate  system  is  defined  by 
the  linear  sensor  array  20  carried  by  the  tow  body  10  and  carries 
the  subscript  "a".  The  origin  is  coincident  with  the  center  of 
the  acoustic  aperture.  The  x-axis  (xa)  is  perpendicular  to  the 


line  passing  through  the  sensors  of  the  linear  sensor  array  20  and 
oriented  in  the  direction  of  motion  when  the  tow  body  10  is 
operated  in  straight  and  level  flight.  Level  flight  in  this 


context  is  that  the  xa_ y a  formed  by  the  coordinates  of  the  linear 


sensor  array  20  is  parallel  in  the  xnv-ynv  plane  of  the  vehicle 
carried  NED  coordinate  system. 
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[0077]  The  y-Axis  (ya)  is  co-linear  with  the  line  passing 


through  the  sensors  22  and  oriented  to  starboard.  The  z-Axis  (za) 

defined  by  a  right  hand  coordinate  system  is  oriented  down  in  flat 
and  level  flight.  The  position  vector  Pa  for  a  location  in  the 


array  system  is  denoted  as 


[0078] 


P  = 


(8) 


[0079]  Projection  of  a  beam  pattern  onto  the  seafloor  500 
begins  with  a  series  of  coordinate  transformations.  In 
particular,  the  set  of  seafloor  locations  that  are  contained 
within  an  area  where  the  linear  sensor  array  20  would  be  capable 
of  detecting  a  ULB  signal  are  transformed  from  the  local  NED 
system  of  coordinates  into  the  array  coordinates  by  the  algorithms 
that  follow.  Following  this  coordinate  transformation;  the  beam 
pattern  is  projected  onto  the  seafloor  500  where  the  result  is 
interpreted  as  a  likelihood  function  used  to  represent  the  area  of 
the  seafloor  from  which  the  detected  ULB  signal  was  likely  to  have 
originated. 

[0080]  Seafloor  grid  locations  are  first  transformed  from 
positions  Pn  in  the  local  NED  system  to  positions  Pnv  in  the 
vehicle  carried  NED  system  of  coordinates  as 
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[0081]  Pnv~Pn~PAn  (9) 

where  PAn  is  the  position  of  the  linear  sensor  array  20  in  the 
local  NED  system. 

[0082]  Seafloor  grid  locations  are  then  rotated  from  their 
positions  in  the  vehicle  carried  NED  system  into  the  array  system 
of  coordinates.  The  array  system  of  coordinates  is  centered  on 
the  acoustic  source  array  and  is  illustrated  in  FIG.  4  and  the 
array  system  is  also  illustrated  in  FIG.  7  with  the  other 
coordinate  systems  used. 

[0083]  The  angular  orientation  of  the  linear  sensor  array  20 
relative  to  the  vehicle  carried  NED  system  is  defined  by  the  set 
of  Euler  angles  that  describe  the  yaw,  pitch  and  roll  of  the  tow 
body  10.  The  Euler  angles  are  provided  by  the  magnetic  heading 
and  gravity  sensors  that  are  included  in  the  tow  body  10.  The 
Euler  angles  change  with  time  and  are  detected  by  the  magnetic 
heading  and  gravity  sensors. 

[0084]  The  yaw  angle  V nv  is  the  angle  from  the  vehicle  carried 
NED  xnv  axis  to  the  projected  vector  of  the  array  xa  axis  on  the 

xnv~Ynv  plane  of  the  vehicle  carried  NED  frame.  The  right  hand 
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rotation  is  taken  about  the  vehicle  carried  NED  znv  axis.  The 
vehicle  carried  NED  frame  is  transferred  to  a  once  rotated 
intermediate  frame  by  a  rotation  matrix  as  =  Pnv^il<r-nv  where 


[0085] 


*/  1< — nv 


cosy/nv  sin* i//nv  0 
-sin  \f/nv  cosy/nv  0 
0  0  1. 


(10) 


[0086]  Each  rotation  is  treated  separately  to  make  following 
the  coordinate  transformation  easier  to  follow. 

[0087]  The  pitch  angle  0nv  is  the  angle  from  the  Xu  axis  of  the 


once  rotated  intermediate  frame  to  the  array  frame  Xa  axis.  The 

right  hand  coordinate  rotation  is  about  the  Yu  axis  of  the  once 
rotated  intermediate  frame.  The  once  rotated  intermediate  frame 
is  transferred  to  a  twice  rotated  intermediate  frame  by  a  rotation 
matrix  "R"  as  =  <  il  where 


[0088] 


Rt2<-i] 


cos0nv  0  -sin 0nv 
0  1  0 
sin 0nv  0  cos 6nv 


(11) 


[0089]  The  roll  angle  Vn  is  the  angle  from  the  Yi2  axis  (or  Zi2 
axis)  of  the  twice  rotated  intermediate  frame  to  that  of  the  body 
frame.  The  twice  rotated  intermediate  frame  is  transferred  to  the 
array  frame  by  a  rotation  matrix  as  Pa  =  ^2^ a <  j'2  where 
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10  0 
0  cos  <pnv  sin  (pnv 

0  —sin  (pnv  cos (pnv 

[0091]  Finally,  the  positions  on  the  seafloor  500  in  the  array 
(Cartesian)  system  of  coordinates  are  converted  to  the  spherical 
system  of  coordinates  illustrated  in  FIG.  4  and  described  in 
Equations  (3)  and  (4)  as 

[0092]  ra  =  +yl  +Zl  ,  ^«=tan  1  and  6a  =  cos  1  ~~  (13) 

\xaJ  \ra 

where  ra  is  the  distance  from  the  origin  of  the  array  coordinate 

system  to  the  seafloor  grid  location  ^ a  \Xa>ya’Za\  with  T 

being  the  transpose  operator.  For  this  invention,  the  transpose 
operator  is  a  standard  operator  in  linear  or  matrix  algebra. 

The  transpose  operator  converts  the  row  vector  ~  [^hz  ’ ^ a  ] 
to  a  column  vector  as  needed  to  complete  the  calculations.  Pa  is 
the  position  vector  that  points  from  a  center  of  the  linear 
sensor  array  20  to  a  location  on  the  seafloor. 

[0093]  A  likelihood  function  for  the  location  of  the  ULB  is 
then  generated  by  projecting  the  detecting  beam  pattern  (s)  onto 
the  seafloor  500.  First,  the  domain  of  the  seafloor  500  from 
which  the  linear  acoustic  sensor  array  20  could  possibly  detect  a 
ULB  is  identified.  In  practice,  it  is  advantageous  to  project  the 
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beam  pattern  onto  the  seafloor  500  at  some  multiple  of  the 
detection  range  estimated  by  evaluation  of  the  passive  sonar 
Equations  (1)  and  (2) .  Suitable  multiples  range  from 
approximately  two  to  four  to  ensure  that  the  beam  projection 
includes  the  location  of  the  ULB,  plus  some  margin. 

[0094]  The  margin  is  based  on  the  estimated  detection  range  for 
the  system  and  will  be  influenced  by  the  local  ambient  noise 
intensity  and  also  by  the  nature  and  arrangement  of  wreckage  on 
the  seafloor.  For  example:  the  ULB  could  be  buried  in  sediment 
and  shielded  by  aircraft  wreckage;  thereby  making  the  ULB  more 
difficult  to  detect.  There  is  no  direct  measure  of  this  margin. 
[0095]  The  identified  seafloor  coordinates  are  then  transformed 


from  the  local  NED  Cartesian  system 


nv> 


to 


the 


spherical  system  of  array  coordinates^,  I/a  ’  ^ a  ’  ^ a  \  by  Equations 
(9)  through  (13) . 

[0096]  The  detecting  beam  pattern  is  then  calculated  at  the 
angular  directions  V^a^a  over  the  set  of  seafloor  grid  points 
that  reside  within  the  determined  multiple  of  the  passive  sonar 
maximum  detection  range.  The  seafloor  projection  of  the  detecting 
beam  pattern  is  then  normalized  such  that  the  sum  over  the 
projected  area  of  seafloor  is  equal  to  one.  The  normalization 
occurs  after  the  procedures  in  the  succeeding  paragraphs  have  been 


28  of  31 


Attorney  Docket  No.  300087 

applied.  Then  the  numerical  values  of  the  beam  pattern  projected 
onto  the  seafloor  are  divided  by  the  sum  of  all  of  the  (non-zero) 
numerical  values  projected  onto  the  seafloor.  This  normalized 
beam  projection  map  is  then  added  to  the  maximum  likelihood 
estimate  which  is  the  sum  of  all  previous  normalized  beam 
projection  maps.  Thus,  the  beam  projection  is  interpreted  as  a 
likelihood  function  that  describes  the  location  of  the  ULB. 

[0097]  The  value  of  the  likelihood  function  for  the  location  of 
origin  of  a  detected  ULB  signal  is  equal  to  either  zero,  or  is 
proportional  to  the  value  of  the  detecting  beam  pattern  at  the 
angular  location  in  array  coordinates.  Also,  the  likelihood 

that  a  detected  ULB  signal  originated  from  a  location  not  on  the 
seafloor  is  equal  to  zero.  Furthermore,  the  likelihood  that  a 
detected  ULB  signal  originated  from  a  location  on  the  seafloor, 
and  at  a  distance  from  the  acoustic  sensor  array  20  that  does  not 
exceed  the  predetermined  multiple  of  the  passive  sonar  detection 
range  is  equal  to  one.  Still  further,  the  likelihood  that  a 
detected  ULB  signal  originated  from  a  location  on  the  seafloor, 
and  at  a  distance  from  the  array  that  exceeds  the  predetermined 
multiple  of  the  passive  sonar  detection  range  is  equal  to  zero. 
[0098]  This  process  is  repeated  with  each  detection  of  a  ULB 
signal,  including  detection  of  the  same  signal  on  more  than  one  of 
the  beams  formed  from  the  array  data.  The  result  is  a  map  of 
seafloor  grid  locations,  provided  in  the  local  NED  system  of 
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coordinates,  representing  the  spatial  likelihood  function  for  the 
location  of  the  ULB. 

[0099]  The  result  of  the  process  is  illustrated  in  FIG.  8, 
where  a  seafloor  grid  502  is  represented  in  kilometers  North  and 
East  of  the  local  NED  system  origin.  The  dotted  line  of  the 
illustration  represents  the  track  taken  by  the  tow  body  10  over 
the  seafloor  500  during  one  leg  of  a  typical  search  pattern. 
Contours  to  the  southeast  of  the  search  track  depict  the 
likelihood  function  for  the  ULB  location.  The  maximum  likelihood 
location  for  the  ULB  is  provided  by  the  "+"  marker  located  at  grid 
position  12.12  kilometers  east  and  7.88  kilometers  north  of  the 
local  NED  origin.  The  likelihood  that  the  ULB  is  located  at 
nearby  locations  diminishes  with  distance  from  the  "+"  marker  as 
shown  by  the  contour  lines  504. 

[0100]  A  noteworthy  feature  of  the  inventive  method  is  that  a 
map  illustrating  an  estimate  for  the  location  for  a  ULB  is 
generated  in  the  course  of  a  single  pass  by  the  tow  body  10.  The 
mapping  uses  data  from  the  linear  sensor  array  20.  It  is  also 
noteworthy  that  spurious  detection  of  potential  interfering 
sources,  such  as  fathometers  or  other  non-stationary  acoustic 
transducers  do  not  accumulate  as  false  targets  because  their 
contributions  to  the  maximum  likelihood  map  are  spatially 
dispersed  in  accordance  with  their  movement. 
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[0101]  The  scope  of  the  present  invention  was  constrained  to 
maintain  compatibility  with  an  existing  system.  Thus,  the  linear 
acoustic  array  20  fits  within  the  geometric  constraints  of  known 
tow  bodies.  The  electronics  required  for  signal  conditioning  and 
communications  will  fit  within  the  volume  provided  by  an  existing 
pressure  vessel  of  the  tow  body  10. 

[0102]  Alternatives  to  present  invention  include  the  use  of 
arrays  with  other  geometries,  such  as  a  planar  array,  but  these 
alternatives  have  greater  complexity  and  expense  compared  to  the 
relatively  short  and  simple  linear  acoustic  array  20  disclosed  by 
this  invention. 

[0103]  It  will  be  understood  that  many  additional  changes  in 
the  details,  materials,  steps  and  arrangement  of  parts,  which  have 
been  herein  described  and  illustrated  in  order  to  explain  the 
nature  of  the  invention,  may  be  made  by  those  skilled  in  the  art 
within  the  principle  and  scope  of  the  invention  as  expressed  in 
the  appended  claims. 
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UNDERWATER  ACOUSTIC  BEACON  LOCATION  SYSTEM 

ABSTRACT  OF  THE  DISCLOSURE 

A  method  using  a  tow  body  and  sensor  array  for  locating  an 
underwater  location  beacon  is  provided.  The  steps  of  the  method 
are  moving  the  body  and  array  in  a  direction  of  motion; 
orientating  the  array  orthogonal  to  the  direction  of  motion; 
beamforming  acoustic  conical  beams  from  the  array  onto  a 
seafloor;  detecting  a  time  series  of  acoustic  data  channels  from 
the  beacon;  converting  the  data  channels  into  spatial  beams  of 
acoustic  frequency;  determining  if  a  signal  of  the  underwater 
location  beacon  is  present  from  each  spatial  beam;  calculating 
likelihood  functions  for  locations  of  the  beacon  to  represent  an 
area  of  the  seafloor  from  which  the  signal  of  the  underwater 
location  beacon  is  to  have  originated;  accumulating  likelihood 
functions  for  the  location  of  the  underwater  location  beacon 
over  the  course  of  a  search  pattern;  and  producing  a  grid  for  a 


beacon  location  based  on  the  likelihood  functions. 
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